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Abstract: The Unaided Single Frequency/Single Epoch Global Navigation Satellite 
System (SF/SE GNSS) model is the most challenging scenario for ambiguity resolution in 
the GNSS attitude determination application. To improve the performance of SF/SE-GNSS 
ambiguity resolution without excessive cost, the Micro-Electro-Mechanical System Inertial 
Measurement Unit (MEMS-JJV1U) is a proper choice for the auxiliary sensor that carries out 
the inertial attitude augmentation. Firstly, based on the SF/SE-GNSS compass model, the 
Inertial Derived Baseline Vector (IDBV) is defined to connect the MEMS-IMU attitude 
measurement with the SF/SE-GNSS ambiguity search space, and the mechanism of inertial 
attitude augmentation is revealed from the perspective of geometry. Then, through the 
quantitative description of model strength by Ambiguity Dilution of Precision (ADOP), 
two ADOPs are specified for the unaided SF/SE-GNSS compass model and its inertial 
attitude augmentation counterparts, respectively, and a sufficient condition is proposed for 
augmenting the SF/SE-GNSS model strength with inertial attitude measurement. Finally, in 
the framework of an integer aperture estimator with fixed failure rate, the performance of 
SF/SE-GNSS ambiguity resolution with inertial attitude augmentation is analyzed when the 
model strength is varying from strong to weak. The simulation results show that, in the 
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SF/SE-GNSS attitude determination application, MEMS-IMU can satisfy the requirements 
of ambiguity resolution with inertial attitude augmentation. 

Keywords: GNSS attitude determination; SF/SE-GNSS compass; inertial attitude 
augmentation; IDBV; MEMS-IMU; ADOP; integer aperture estimation 



1. Introduction 

The weak Global Navigation Satellite System (GNSS) model is defined as a scenario with 
insufficient observation quantity or poor measurement quality, so that the availability of the GNSS 
application becomes worse. As to the GNSS attitude determination, weak GNSS model strength 
can be caused by few available GNSS satellites and frequencies, insufficient epochs and noisy 
measurements, etc. [1-5]. These factors lead to the failure of reliable but unaided ambiguity resolution. 
The increasing number of urban canyon environments makes satellites' signals easily blocked or 
contaminated by multipath signals. The highly dynamic movement of vehicles also indirectly weakens 
the model strength [6]. In these situations, it is difficult to realize the ambiguity resolution of the GNSS 
attitude determination system based on low-end receivers. 

Unaided SF/SE-GNSS ambiguity resolution is a typical GNSS application in a weak GNSS model. 
The augmentation methods can be categorized into three classes: multi-GNSS, multi-sensor 
information fusion and a priori constraints. Multi-GNSS can increase the data redundancy and 
improve the system reliability [7]. Multi-sensor information fusion, or integrated navigation, can 
overcome the shortcomings of different sensors. The vehicle attitude, baseline vector or difference 
carrier phase measurements are usually input to the optimal estimator for improving the availability of 
GNSS applications with weak model strength [8,9]. Typically, a priori constraints include the baseline 
length [10], the geometry of multi-baselines [11] and the non-holonomic constraints [12]. All these 
constraints can be used to decrease the ambiguity search space and improve the success rate of 
ambiguity resolution. 

From the 1990s on, methods of inertial measurement augmentation for ambiguity resolution 
appeared successively in both studies on GNSS attitude determination and Strapdown Inertial 
Navigation System/GNSS (SINS/GNSS) integrated attitude determination. In 1998, to improve the 
accuracy of GNSS attitude determination, Hayward and Egziabher made use of gyro measurements to 
smooth the high frequency errors contained in GNSS carrier phase measurements [13,14]. In 1999, 
Han showed that the inertial attitude measurement from a rate-gyro can be used to improve the success 
rate of ambiguity resolution [15]. From 2001 to 2003, Yang revealed that the MEMS-IMU 
measurement can be used as an a priori constraint to improve the integrity of ambiguity 
resolution [9,16]. In 2004, to resolve the problems existing in the GNSS attitude determination 
application when low-end receivers are used, Wang introduced the gyro-rate measurement into 
the attitude integration filtering and improved the performance of GNSS attitude determination 
system [17,18]. In 2008, Dickman studied how to use the inertial measurement to depress the low 
frequency errors contained in GNSS carrier phase measurements [19], but this method requires a high 
IMU measurement precision, hence, it is not applicable to GNSS attitude determination. 
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In theory, the method of inertial attitude augmentation for ambiguity resolution is not the same as 
that of the multi-sensor information fusion augmentation. Traditional SINS/GNSS integrated 
navigation systems use the SINS measurement error as the state vector, and the GNSS output as the 
observations. However, the error characteristics of low-end IMU measurements vary greatly with the 
vehicle dynamics and working temperature. In a weak GNSS model, it is very difficult to achieve the 
optimal estimation for SINS errors and the filter divergence will lead to the unavailability of the GNSS 
application. As the output frequency of the low-end receivers gradually increases, the GNSS attitude 
determination system can satisfy the frequency requirements in common dynamic situations. The 
MEMS-IMU measurement within a short span can be used as a high quality prior constraint to 
augment the GNSS model strength. 

In this paper, Section 1 reviewed the development of inertial augmentation for the GNSS model and 
introduces the meanings of this research. Section 2 will give the definition of IDBV and generate an 
inertial augmentation ambiguity search space based on the IDBV and SF/SE-GNSS compass. Then, 
the geometry characteristics of this search space are also analyzed. Section 3 will reveal the 
mechanism of inertial attitude augmentation for SF/SE-GNSS model strength. The ADOP is 
introduced and utilized to analyze the effect of inertial attitude measurement for SF/SE-GNSS model 
strength augmentation. Section 4 will analyze and compare the performance of ambiguity resolution 
based on an integer aperture estimator with fixed failure rate after simulation experiments. Finally, 
Section 5 will provide a summary of the work. 

2. Inertial Attitude Aided Ambiguity Search Space 

This section gives the definition of IDBV, and then generates an ambiguity search space based on a 
float ambiguity vector and its Variance-Covariance Matrix (VCM). The obtained space is called the 
inertial attitude aided ambiguity search space. By analyzing the geometry characteristics of this space, 
how the inertial attitude measurement imposes its influence to the unaided SF/SE-GNSS ambiguity 
search space will be clarified. 

2.1. IDBV 

In a certain GNSS epoch, the attitude measurement of SINS is denoted as C n b , thus, the IDBV can 
be defined as: 

D:=C b V (1) 

In Equation (1), b b denotes the true value of the baseline vector in the body-frame (b-frame). Then, 
£>(•) is introduced to represent the deviation between the measurement and true value. The following 
equation is thus obtained [20]: 

5b n =-(y n x).b n (2) 

In Equation (2), y n =[V",7°,7°] is to represent the rotational angle error of C£ in navigation 

frame (n-frame). Equation (2) establishes the relation between the measurement errors of inertial 
attitude and IDBV. 
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The differential equation of y n with respect to time can be expressed as [21]: 

7 n =-C b "5<-(<x)y" + 5< (3) 

In Equation (3), 8e> b b is the measurement error vector of the gyros, and if MEMS gyros with the 
current accuracy level are utilized, each component of 5co b b has an order of magnitude of 
30~ 200deg/h («1.5xl0 -1 ~lmrad/s). to? denotes the true value of the n-frame angular velocity 
vector with respect to the inertial frame (i-frame) and represented in n-frame. For a common land 
vehicle maneuvering on the Earth's surface, e.g., the maximum linear velocity is limited to 200 km/h, 
each component of co" n will have the same order of magnitude as the angular rate of the Earth's 
rotation, i.e., Q. = 7.2921 15 xlO 2 mrad/s . After a cross product with y n , each component of which is 
assumed to be less than 0.01 deg ( « 1.7x10 'mrad), the second item on the right side of Equation (3) 
will have a 1.2x10 2 mrad/s order of magnitude, dco^ denotes the measurement error contained in the 
measurement of oi° . If the time intervals between different GNSS observations are very short, and 
GNSS position and velocity are both available with accuracy levels of 10 m and 0.1 m/s, respectively, 
and it can be proved that each component of 5(0^ has a 10 4 mrad/s order of magnitude. Hence, 
Equation (3) can be approximated as: 

r - -c»M b (4 ) 

It can be noted that Equation (4) omits the influence of vehicle dynamic. Hence, in this situation, 
the measurement error of the gyros is the principal contributor to the divergence speed of the inertial 
attitude measurement error. 

We denote the neighboring two epochs as t 0 and t x , t 0 is the starting time of the SINS calculation. 
The initial attitude at t 0 is provided by a GNSS attitude determination system with two orthogonal but 
equal-length baselines. The inertial attitude will be the integral of the MEMS-IMU measurements and 
starting from the initial attitude. With Equation (4), the inertial attitude measurement error at t y equals: 

y n (O=Y n (O-JX n (0 8 <(0^ (5) 

The common model of Scoj^ is [21,22]: 

K b =K b Bia ;+w b Ran ; (6) 

where K b ias ' is the equivalent gyro bias, which is usually modeled as a random constant vector. The 
equivalent measurement noise w Rand ' is modeled as a white noise vector. Substituting Equation (6) into 
Equation (5), then it is obtained that: 

r (0 = r ) - £ c a b {t)[Kj ) dt - £ c b - (t) ( w b Ran ; ) * 

= y n (f 0 ) - £ (0 ( K Bias' ) dt - £ ( W" Rand ' ) dt (7) 
= y n ^)-£Q(0( K B ia ;)^-v n Walk 

C b is an orthogonal matrix, hence, w Rand ' = C b (V)^w Rand ' j retains the statistical characteristics of 
w Lnd [12,23], which means w Rand ' is still a white noise vector. A random walking item v^, alk will be 
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obtained by integrating the white noise vector w^ and ' with respect to t . On the other hand, C£ (f) is 
continuous in closed interval [? 0 ,^] and integrable in open interval (t Q ,t l ) . According to the first 
integral mean theorem, there is: 

£q(0(kL')^-c:(^)£(kL')^ 

-c:(^)(K b Bia ;)(^-^)-(KB ia ;)^ 

In Equation (8), there is £e(f 0 ,f,) and K° ias ' in Equation (8) retains the stochastic characteristics 



(8) 



of K 



b ' 

Bias ' 



According to Equations (2), (7) and (8) and the error propagation rule, the VCMs of y n and b" can 
now be derived: 



Q „ =cr!»I 3 +Q 



Ya, 



Yo ' 



<4 =^Bias^ 2 +< nd Ar, 



(9) 



Q/=cov(K) = (b n x)Q y „ (b n x) T 



where At = t l -t 0 , <j^ ias and cr^ and denote the variances of the components contained in Kg ias and 
w Lnd • Q ? » 1S tne VCM of the initial attitude measurement provided by the GNSS attitude 

determination system. Assuming that the GNSS attitude determination system has the equal length but 
orthogonal dual-baseline style, then Q ,, can be written as [20]: 



^YS b 2 





"0" 




<2i 3 -c; 


0 


[o o i]c*. 




1 





(10) 



where b is the baseline length, a\ denotes the variance of the GNSS carrier phase measurement. 
Substituting Equation (10) into Equation (9), then we obtain: 



0 



0 



0 bW +al 

Jtt 

0 0 



0 
0 



b 2 a\ +2ai 



it, 



(11) 



Equation (11) is derived basing on the presumption of b b = [Z?,0,0] r . For convenience in the 
following analysis, Q/ can be enlarged to: 

Q l =(b 2 ^ yl +2al)l 3 =afl 3 (12) 
Comparing with Q t ' , Q : evaluates the precision of b" more conservatively. 
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2.2. The Geometry Characteristics of the Inertial Attitude Aided Ambiguity Search Space 

The GNSS attitude determination system with low-end receivers usually makes use of the double 
difference carrier phase measurement model [20], its SF/SE-GNSS counterpart is [24]: 

y=iD r -A€» = D r Ab + ia + n, bel 3 , aeZ" 4 (13) 

where the dimensions of y and AO are m-1 and m respectively; D r is a double difference operator 
matrix with the dimension (m - 1) x m; A is constructed by the Line-Of-Sight (LOS) unit vectors 
between user and the available satellites; b denotes the baseline vector. 
The VCM of y is given as: 

Q 7 =o&m. <rL=2*l (14) 

The standard deviation of un-differential GNSS carrier phase observation cr^ has the millimeter 
order of magnitude [25,26]. 

Substitute b" into Equation (13) and a float solution vector a can be resolved: 



a^^y-D^Ab; 1 ), 



(15) 



Q 4l =-^(Q y +D r AQ I A r D). 

Jlj is the so-called float solution of Inertial Derived Ambiguity Vector (IDAV). 
With a t and , a double difference ambiguity search space can be written as: 

T I (a) = (a I -a) r Q- 1 (a I -a)</ (16) 

The ambiguity search space described by Equation (16) is a (hyper)ellipsoid. Both its center and 
shape vary with the reference satellite selection. Therefore, before analyzing the geometry 
characteristics of Tj(a)<;}f 2 , it is necessary to lift the standard double difference grid aeZ m_1 into 
the single difference ambiguity search space jg- [27], then an ambiguity search space that is 
independent of the selection of reference satellite can be obtained: 

T^sXsl-sfDQl^fsl-s')^ 2 , 

(17) 

a,=D r s;, a = DV. 

As the equivalent form of Equations (16) and (17) can be regarded as the standard form of the 
inertial attitude aided ambiguity search space. 

The shape of the (hyper)ellipsoid T t (s') < % 2 is mainly determined by DQ s 'D r . Referring to the 

proof of theorem 1 in [27], if m is larger than 4, the eigenvalue decomposition of DQ^'D r is given as: 

DQ-D r =RAR r , 



A = diag 
R = [w V U]. 

where: 



V 2cr o> 2cr o 



(18) 
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T = diag[y l (s l ) y 2 (s 2 ) r 3 (s 3 )], 

w = e ffl (e>,„)" /2 , V=E(E r E)" /2 , (19) 
U = F(F r F)' /2 , F = PAS r , S = [s, s 2 s 3 ]. 

In Equation (19), e m is the m-dimension vector with all elements equal to 1; E is a basis matrix of 
R(D)DN(A r P), where N(A r P) describes the null space of A r P and R(D) means the column 
space of D; P = D(D r D) D r is the orthogonal projection operator of R(D) ; The baseline gain 
y i (s,. ) , i = 1, 2, 3 is defined as: 

Different from the baseline gain number defined in [24], the gain number defined by 
Equation (20) describes that in the direction of gain vector s ; eK 3 , the accuracy of the baseline 
solution b(O) varies with respect to b" . Q 6 (®) is the VCM of b(O). It can be proved that if PA is 
full column rank, s, will be the unit eigenvector for the eigenvalue decomposition of 
A r PA = (PA) r ■ (PA) , which means: 

A r PA = SES r , 

S = [s, s 2 s 3 ], (21) 

Since A r PA is positive definitive, the three eigenvalues in Equation (21) satisfy that 

^. >0, i =1,2,3,. Then, it can be proved that the baseline gain /,-(s,-), i = 1,2,3 are the functions of 
£>0, i=l,2,3: 

YM,) = ^T^ i = W (22) 

Hence, with Equation (17) in Equation (22), the decoupled expansion of the inertial attitude aided 
ambiguity search space can be written as: 



T, 



i^«-oT*S[<«-')T}^ (23, 



From Equation (23), it is obvious that the principal axis lengths in the V-axis of the (hyper)ellipsoid 
Tj (s') < % 2 are all equal to -Jl^a^/A, and the three principal axis lengths in U-axis are equal to: 



lu i =J r ~+i- X °* =^ L ^Ratio_IG 2 -% i + 2, i = 1,2,3 



X 



(24) 



In Equation (24), Ratio _IG-a l /a <S) ; ,i = 1,2,3 are not only the eigenvalues of A r PA , but also 
the singular values of PA . PA can be seen as the projection of the user-satellite geometry matrix A 
onto R(D) . Thus, <^-,z = 1,2,3 describe the geometry of the used satellite. 
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3. Analysis for Augmenting the SF/SE-GNSS Model Strength with Inertial Attitude Aided 

From the double-difference measurement model of the SF/SE-GNSS compass, the ambiguity float 
solution vector and its VCM can be resolved: 

a = |[y-D r Ab(0,p)", 

i " (25) 

Q fi =-[2^D r D + D r AQ fi (p)A r D] 

In Equation (25), both b(0,p) and a are the real-valued estimations of the double-difference 
measurement model of the SF/SE-GNSS compass, where (®,p) shows that the baseline vector 
estimation b uses both code and phase observations. The necessary condition for Equation (25) is the 
unique solution is that the number of available satellites is more than 4. The standard form of the 
ambiguity search space determined by a and Q- is: 

T(s') = (§' -s' f DQ- D r (§' - s') < z 2 (26) 

Comparing Equation (17) and Equation (26), it can be proved that DQ-'D r has the same 
eigenvectors as those of DQ-W, i.e., w, V and U ; the eigenvalues corresponding to w and V are 
equal to those in the DQ-^D 7 case, respectively; the eigenvalues of U = (UpU 2 ,u 3 ) are equal to 
X 1 j2(cT 2 p + <r 2 ) [27], or the axis lengths in the directions of U -axis of (hyper)-ellipsoid T(s') < j 2 are 
all equal to: 

lu = ^ Ratio _ PC 2 -2 + 2 (27) 

where Ratio _ PC = <j p l<j< s> , a p is the standard deviation of code observation. 

Comparing Equation (27) with Equation (24), the variations caused by inertial attitude measurement 
will be reflected in the directions of U -axis, and lu is independent of the user-satellite geometry, 
which is revealed by <^ =2 in Equation (27). If Ratio _IG 2 ■ £ )i < Ratio _ PC 2 -2 , or 

a 2 <a 2 -(2/^.),/ = 1,2, 3 , there will be lu^lu' , thus it means that the volume of T t (s')<j 2 is 

smaller than that of T(s') < j 2 , or the inertial attitude measurement decreases the unaided SF/SE-GNSS 

ambiguity search space. 

The magnitude of can be revealed by the simulation results. Firstly, based on the GNSS Yuma 

file, the positions of GPS and BD2 satellites within 12 h can be calculated. Then, setting the user's 
position on Earth, the maximum number of available satellites can be determined, and the geometry 
matrix A is constructed. Finally, = 1,2,3 can be obtained by the eigenvalue decomposition of 

A T PA at each epoch. In this simulation, the user position is set as one point in Changsha, 28.2202°N, 
112.9925°E, and the number of available satellites and A are recorded every 6 min. At each recording 
epoch k , the maximum eigenvalue of A r PA is denoted as ^ -1,2,3 . The variation of £ max 

during 12 h and the mean values are shown in Figure la. 
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Figure 1. (a) The variation of £ max and the mean values during 12 h in Changsha; 
(b) The variation of maximum numbers of available satellites in Changsha during 12 h. 

7 1 1 1 i — 



s> 



o GPS max values 

— GPS mean values 
G BD2 max values 

— BD2 mean values 




16 



4 6 
Time/hour 

(a) 



2 
w 

o 

n 

a 

> 

CO 



mmuj i.i.i.i.i.t i i.i I.U.L _t. 

H — I ff wSsSU hm9W w 



E 



11 
10 



.- # 
x 

CO 

E 
o 



* GPS 
■* BD2 



pi 



■ ■ ■ 



* I.UAI.I.I.I. U.I.I.I.I. 1.1.1.1 1.1.1.1.1.1 

^tJAOJ^. jJjJAfcf... ^-j-j-j-j-j-j-lif 
■Itil-iiiiiiitl-- "lititii" JJik. ^ILii-i-fr 

^T 7 -TT 



10 



12 



Time/hour 

(b) 

It must be noted that to calculate <f max at each recording epoch, all the available satellites are used. 
Hence, the corresponding maximum numbers of available satellites are shown in Figure lb. Above all, 
for the inertial attitude measurement decreasing the SF/SE-GNSS ambiguity search space, a sufficient 
condition is that: 



2 

a-< 



maX (Cax) P 



(28) 
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Equation (28) can be used as an empirical formula for selecting the inertial sensors for inertial 
attitude augmentation. Until now, the analytical relation between the inertial attitude measurement and 
the SF/SE-GNSS ambiguity search space has been built, and the influences of the inertial attitude 
measurement to the unaided SF/SE-GNSS ambiguity search space are directly revealed. From the 
perspective of geometry, the mechanism of inertial attitude augmentation is revealed for the first time. 

To describe the model strength of SF/SE-GNSS ambiguity resolution, the ADOP is introduced. This 
subsection will quantitatively describe the augmentation effect based on inertial attitude measurement 
with different ADOPs. The formula defining ADOP is given by [28]: 



ADOP = l ' m ^det(Q k ) 



(29) 



where m - 1 is the dimension of Q- . ADOP has the same unit as ambiguity, i.e., cycle. 

According to Equation (29), the determinant of VCM of the ambiguity float solution vector is 
needed to calculate ADOP. Based on Equation (15), the determinant of Q- equals to: 



2a, 

X 2 



2 \ 



D r D + — -D r AQ,A r D 

2al 



Utilizing the determinant decomposition formula [29]: 

|e - fg h| |g| = |e| |g - he _1 f| 

Together with Equation (31), Equation (30) can be rewritten as: 



QJ = 



2a\ 

v ^ , 



|D r D||l 3+ 0^(0)1 



(30) 



(31) 



(32) 



In Equation (32), it is easy to obtain that |D r D| = m. According to Equations (12)— (14), it can be 
proved that QjQ: 1 (O) = SrS r , thus: 



m 



/ t \ m-1 

' 2a ^ 3 



f[ U + jRatio_IG 



(33) 



The same process can be imposed on Equation (25). Then the determinant of Q a equals: 



\QA = m- 



2a 



2 \ 



m-1 



1 ^ 

l+ — Ratio PC-2 
2 



(34) 



Substituting Equations (33) and (34) into Equation (29) 



ADOP, 



2a^ 

X 2 



Yl\ 1 + -• Ratio _IG 2 ^ 



(35) 
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In Equation (35), ADOP l and ADOP SE/SE represent the ADOPs of the inertial attitude augmentation 
model and unaided SF/SE-GNSS model, respectively. According to Equation (35), it is known that, if 
a certain accuracy level of GNSS measurement was given, ADOP SF/SE would only be determined by 
the number of available satellites, and the user-satellite geometry has no effect on ADOP SF/SE . Besides 
this, it also can be derived that Equation (28) is a sufficient condition for augmenting the SF/SE-GNSS 
model strength with the aided inertial attitude. 

Considering that the number of available satellites equals four, five and six, the simulation results 
indicating how the ADOP SF/SE values vary with respect to cr^ and Ratio _ PC , are explicitly shown in 
Figure 2. 

Figure 2. The variation of ADOP SF/SE when the number of available satellites equals four, 
five and six, respectively. 
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In Figure 2, the varying ranges of and Ratio _ PC are both determined by the actual performance 
of several GNSS receivers, which are listed in Tables 1-3. 



Table 1. The nominal measurement precisions of three GNSS OEM cards. 





Frequency 


Code(lcr) 


Carrier Phase ( \a ) 


Ratio _ PC 


Novatel OEM628 


GPS-L1 


4 cm 


0.5 mm 


80 


Novatel OEMStar 


GPS-L1 


5 cm 


0.6 mm 


83.3 


UB240 


GPS-L1/BD2-B1 


10 cm 


0.5 mm 


200 
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Table 2. Standard deviation of BD2 code and phase observations in static experiments. 



C07 
(Reference SV) 


Cll 


C14 


C04 


C01 




(mm) 


N/A 


1.2864 


1.0368 


1.1712 


0.8448 




G p (cm) 


N/A 


20.04 


13.3 


16.29 


11.21 




Ratio _ PC 


N/A 


155.8 


128.3 


139.1 


132.7 






C03 


C06 


C02 


C05 


C09 




V® (mm) 


0.8064 


0.7872 


1.1904 i 


0.8064 


1.0176 




G p (cm) 


10.78 


10.55 


14.64 


10.66 


12.77 




Ratio _ PC 


133.7 


134 


123 


132.2 


125.5 




Table 3. Standard deviation of GPS code and phase observations in static experiments. 


G24 G15 
(Reference SV) 


G18 


G14 




G22 


G21 


(mm) n/A 


0.8373 


0.5899 


1.2369 


1.1037 


0.7231 


o p (cm) n/A 


9.67 


9.5 


18.66 




13.53 


9.68 


7?a//o _ PC n/A 


115.5 


161 


150.9 




122.6 


133.9 



This paper will implement the ambiguity resolution basing on the integer aperture estimator with 
fixed-failure rate [30,31]- This estimator requires that the success rate of ambiguity resolution be larger 
than 0.99, which corresponds to the interval ADOP< 0.15. Hence, the reference level of ADOP 
appeared in Figure 2 was set as 0.15. 

The simulation results shown in the first (upper) subplot of Figure 2 indicate that when 
ADOP SP/SE < 0.15 , while there are four available satellites, the values of Ratio _ PC cannot be 
achieved by current GNSS receivers. If the number of available satellites is five and >2 mm, i.e., 
the green line and the yellow line shown in the second (bottom, left) subplot of Figure 2, the values of 
Ratio _ PC make ADOP SP/SE <0.15 are still hard to satisfy with current low-end GNSS receivers. As 
to six available satellites and (7^ < 1 mm, the ADOP SFISE values in the given interval of Ratio _ PC 
will always be less than 0.15, just like the results shown by the red line in the third (bottom, right) 
subplot of Figure 2. 

To weaken the model strength of the SF/SE-GNSS compass with aided inertial attitude 
measurement, in the following simulation, is set to be 3 mm and let be equal to 3. The precision 
combinations for generating the MEMS gyro measurements are shown in Table 4. 

Table 4. The simulation parameters of MEMS gyro. 



Gyro Bias ( \a ) Gyro Noise( \a ) 



Combination 1 


30°/h 


10%/h 


Combination 2 


80°/h 


30°/^ 


Combination 3 


100°/h 


5()7Vh 


Combination 4 


200°/h 


80°/^ 
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The MEMS accelerator bias and noise level are chosen as 5x10 3 g/h and 1x10 3 g/h , respectively. 
Since only the extra-short baseline scenario is considered, the baseline length is set to be successively 
0.5, 1.0, 1.5 and 2.0 m. 

When the number of available satellites successively equals four, five and six (Figure 3), the 
simulation results for ADOP l are generated under different conditions of baseline lengths and MEMS 
gyro precision combinations. 

Figure 3. (a) The variation of ADOP l with four satellites available; (b) The variation of 
ADOP l with five satellites available; (c) The variation of ADOP l with six satellites available. 



Combination 1 



Combination 2 



0.15 




0.15 



0.15 



36 47 70 

GNSS Epoch/Second 



Combination 1 



(a) 



19 24 36 70 90 

GNSS Epoch/Second 



Combination 2 




291 387 580 

GNSS Epoch/Second 




0.15 



109 146 218 

GNSS Epoch/Second 
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Figure 3. Cont. 



Combination 3 



Combination 4 




0.15 



117 175 
GNSS Epoch/Second 



175 

GNSS Epoch/Second 



(b) 



Combination 1 



Combination 2 




729 971 1456 

GNSS Epoch/Second 

Combination 3 



O 
□ 
< 



0.15 




220 292 438 

GNSS Epoch/Second 



O 
□ 
< 



0.15 




274 365 547 

GNSS Epoch/Second 

Combination 4 



O 
a 
< 




111 147 220 

GNSS Epoch/Second 



(c) 



In Figure 3a, where only four satellites are available, if the worst accuracy combination 4 is chosen, 
the state of ADOP l < 0.15 can only be kept about 1 min. Otherwise, if the best accuracy combination 1 
is used, the time span for ADOP l <0.15 is still no more than 2 min. Hence, Figure 3a demonstrates 
that four available satellites will be a very weak condition for the model strength of SF/SE-GNSS 
compass with aided inertial attitude measurement. In Figure 3b, where five satellites are available and 
the baseline length equals to 1.5 or 2.0 m, the worst accuracy combination 4, or the second best 
accuracy combination 2 still cannot keep ADOP l <0.15 for more than 3 min. In Figure 3c, whatever 
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the accuracy combination is chosen, the MEMS-IMU attitude augmentation can already keep 
ADOP l < 0.15 more than 3 min, meaning the SF/SE-GNSS model with aided inertial attitude has 
some strength when six available satellites are considered. 

Combining the simulation results shown in Figures 2 and 3, the number of available satellites is the 
principal factor that influences the SF/SE-GNSS model strength. Based on the current accuracy levels 
of the low-end GNSS receivers and MEMS gyros, the key to realize the inertial attitude augmentation 
for SF/SE-GNSS ambiguity resolution is to resolve the existing problems when the number of 
available satellites is between four and six. 

4. Performance Analysis for the SF/SE-GNSS Ambiguity Resolution with Aided Inertial Attitude 

In this section, a specific simulation experiment will be designed with the low-end GNSS receivers 
and the MEMS-IMU considered. In the framework of an integer aperture estimator with fixed-failure 
rate, the performance of SF/SE-GNSS ambiguity with aided inertial attitude will be analyzed based on 
the simulation results, while the rate ratio test and the difference test will be compared under different 
model strength scenarios. 



4.1. Simulation Experiment Design 

Figure 4 explicitly describes the basic simulation flow chart. 

Figure 4. Simulation flow chart of inertial attitude augmentation for SF/SE-GNSS 
ambiguity resolution. 
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The main functions of each module in Figure 4 are introduced here: 

(1) Module A. Assume that the geometric center of primary antenna A and the origin of b-frame 
overlap each other, the true baseline vector in b-frame and initial vehicle attitude are known in 
prior, which are labeled as b b and C£ (? 0 ) respectively. 
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(2) Module B. As to the GNSS data collected from an field experiment, the GNSS position and 
velocity measurements at certain epoch are used as the true initial navigation states, which 
include the time, position and velocity for INS calculation. These parameters are denoted as ? 0 , 
pl(?o) an d v^(? 0 ), successively. 

(3) Module C. The trajectory generator outputs the true values of vehicle position, velocity, attitude 
and IMU measurements, which are denoted as p^ (?) , (?) , C b (?) , f ib b and (o b , and there 
is t > ? 0 . 

(4) Module D. The true satellite positions are calculated from the broadcast ephemeris, which is 
also collected in the same field experiment of Module B, and the true satellite positions are 
labeled as % Yi (?) ,i = 1, 2...m . 

(5) Module E. With the true navigation states utilized, the true baseline vector, satellite positions, 
un-differenced code and carrier phase observations are generated, which are labeled and 
included in y (?) . The true ambiguities are included in the vector a , and the user- satellite unit 
LOS vectors form the design matrix A(?) . 

(6) Module F. After adding the simulation noises into the original outputs of Module C and Module 
E respectively, the MEMS-IMU and GNSS measurements can be obtained and denoted as f£ , 
e> b b and y(?) . According to the standard deviations of the added noises and the setting baseline 
length, the VCM of IDBV can be calculated. 

(7) Module G. Using the MEMS-IMU measurements f b and co b b , the SINS calculation starts 
from the initial navigation states. Besides the SINS navigation states, this module also outputs 
the IDBV b" (?) . 

(8) Module H. Basing on the SF/SE-GNSS float ambiguity vectors given in Equations (15) and (25), 
the integer aperture estimator with fixed failure rate tries to fix the ambiguities. 

(9) Module I. This module compares the ambiguity solutions output from module H with the true 
ambiguity a , and evaluates the performance of SF/SE-GNSS ambiguity resolution with inertial 
attitude measurements aided. 

It should be noted that the influence of the user position error on A(?) can be omitted [32], hence in 
the transmission process from module E to H, module F doesn't add any measurement error to A(?) . 
Considering the double difference operation to 0^ Vi (?) in module H, module F does not model the 
propagation errors and clock errors contained in the original GNSS observation vector y(?) . Hence, 
the noise item n y can be modeled as a white noise vector. Otherwise, in module F, the modeling 
methods of MEMS-IMU angular velocity measurement error are the same as those for 5e) b b in 
Equation (6). 

The integer aperture estimator with fixed failure rate can control its failure rate by setting the testing 
threshold. However, there is no analytical relation between the testing threshold and the failure rate. 
Herein, the threshold will be determined by Monte Carlo simulation. Figure 5 gives the detailed steps 
for the determination of testing threshold in the framework of integer aperture estimator with fixed 
failure rate. In Figure 5, N f denotes the number of failing samples in ambiguity resolution. The 
number of simulations is set as 100,000. 
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Figure 5. Determination of the testing threshold in the framework of integer aperture 
estimator with fixed failure rate. 
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4.2. Performance Analysis 

Since the approximation Equation (4) omits the influence of the vehicle dynamics, the simulation in 
this subsection will only consider the scenario of a static vehicle and be carried out three times based 
on the number of available satellites. The performance analysis for SF/SE-GNSS ambiguity resolution 
also focuses on the influence of model strength varying from weak to strong, i.e., ADOP l ranges from 
0.1 to 0.5. In the framework of an integer aperture estimator with fixed failure rate, the ratio test and 
difference test will be implemented separately. To evaluate the simulation results, the ambiguity 
solutions are classified and summarized in Table 5. 



According to the simulation flow chart shown in Figure 4, each simulation experiment period is set 
as 1000 s. The baseline length is set as 2 m. The accuracy combination 2 shown in Table 4 is chosen 
for the MEMS gyros, and let cr 0 equal to 0.003 m. Based on the broadcast ephemera collected in the 
field experiment, the GPS or BD2 constellation is generated, and the numbers of available satellites for 
GPS and BD2 is equal to eight and 11, respectively. Under any one of the conditions that the number 
of available satellites is limited to be four, five or six, the satellite combination which gives the 
weakest model strength, or in other words, the maximum ADOP l value, is picked out from all the 



Table 5. Ambiguity resolution results classification. 



Style 1 
Style 2 
Style 3 
Style 4 



Accept correct integer ambiguity 
Accept wrong integer ambiguity 
Reject correct integer ambiguity 
Reject wrong integer ambiguity 
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available satellites. The fixed failure rate for integer aperture estimator is set as 0.01 and 0.001. For 
different simulation scenarios, the SF/SE-GNSS ambiguity resolution results with MEMS-IMU 
attitude aided will be given in Figure 6. 

Figure 6. (a) The MEMS-IMU attitude aiding SF/SE-GNSS ambiguity resolution results 
with four satellites available; (b) The MEMS-IMU attitude aiding SF/SE-GNSS ambiguity 
resolution results with five satellites available; (c) The MEMS-IMU attitude aiding 
SF/SE-GNSS ambiguity resolution results with six satellites available. 
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From Figure 6, the output results of integer aperture estimator varied from 1 to 3 or 4, Style 2 never 
appeared in the process for model strength varying from strong to weak. It reveals that when the 
MEMS-IMU attitude measurements are used to augment the SF/SE-GNSS model strength, the integer 
aperture estimator with fixed failure rate has good reliability once the testing threshold is determined. 
Even if the requirement for failure rate was relaxed, the performance of the SF/SE-GNSS ambiguity 
resolution still held at a high level. On the other hand, no matter which testing method was selected 
between the ratio test and the difference test, there was no significant distinction in performance, this 
means that under the condition of SF/SE-GNSS model with aided inertial attitude, the testing method 
is not the primary factor that affects the performance of ambiguity resolution. 

5. Conclusions 

This paper was set up based on the GNSS compass model and the MEMS-IMU attitude 
measurement application. It investigated how to utilize the inertial attitude measurement to enhance 
the performance of unaided SF/SE-GNSS ambiguity resolution. According to the definition of IDBV, 
the relationship between the inertial attitude measurement and the ambiguity search space was 
established. Then, the essence of inertial attitude augmenting the SF/SE-GNSS model strength was 
revealed from the geometrical perspective, while the empirical formula of selecting the inertial sensors 
for inertial attitude augmenting was given. ADOP was introduced to quantitatively describe the model 
strength. Based on the simulation results, the factors influencing SF/SE-GNSS model strength were 
analyzed overall. It was concluded that among the inertial sensor measurement accuracy, the GNSS 
measurement accuracy, the baseline length and the number of available satellites, the last factor is the 
most significant one in practice. In the designed simulation experiment for the SF/SE-GNSS ambiguity 
resolution with MEMS-IMU aided attitude, the integer aperture estimator with fixed failure rate was 
used. During the model strength varying from strong to weak, the performance of MEMS-IMU attitude 
augmenting SF/SE-GNSS ambiguity resolution was analyzed, and some practical suggestions about 
the application of the integer aperture estimator were given. 
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